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We explore the chiral aspects of extrapolation of observables calculated within lat- 
tice QCD, using the nucleon magnetic moments as an example. Our analysis shows 
that the biggest effects of chiral dynamics occur for quark masses corresponding 
to a pion mass below 600 MeV. In this limited range chiral perturbation theory is 
not rapidly convergent, but we can develop some understanding of the behaviour 
through chiral quark models. This model dependent analysis leads us to a simple 
Pade approximant which builds in both the limits m^: — ► and — ► oo correctly 
and permits a consistent, model independent extrapolation to the physical pion 
mass which should be extremely reliable. 



1 Introduction 

At present the only known way to calculate the properties of QCD directly is 
through the formulation of lattice gauge theory on a discrete space-time lattice. 
While the lattice formulation of QCD is well established H, there have recently 
been a number of exciting advances in lattice action improvement which are 
greatly facilitating the reduction of systematic uncertainties associated with 
the finite lattice volume and the finite lattice spacing. However, direct simu- 
lation of QCD for light current quark masses, near the chiral limit, remains 
computationally intensive. In particular, present lattice calculations of masses 
in full QCD are limited to equivalent pion masses of order 500 MeV or more - 
although there have been recent preliminary-results from CP-PACS using im- 
proved actions at a mass as low as 300 MeVB. For nucleon magnetic moments 
the situation is much worse, with the best calculations still being made within 
quenched QCD for pion masses greater than 600 MeV. This is clearly beyond 
the range of validity of chiral perturbation theory (xPT). 



In view of this situation, the present approach of calculating the properties 
of QCD using quark masses away from the chiral regime and extrapolating to 
the physical world is likely to persist for the foreseeable future. It is therefore 
vital to understand the quark mass dependence of hadronic observables calcu- 
lated on the lattice and how to connect these calculations with the physical 
world. A major difficulty in this endeavour is the rapid rise of the pseudoscalar 
mass for small increases in the quark mass away from the chiral limit - with 

oc m, where to is the light quark mass. 

Historically, lattice results were often linearly extrapolated with respect 
to mj, particularly in exploratory calculations. More recently the focus has 
turned to chiral perturbation theory (xPT), which provides predictions for the 
leading nonanalytic quarksmass dependence of observables in terms of phe- 
nomenological parametersQU. Indeed, it is now relatively standard to use xPT 
in the extrapolation of lattice simulation data for hadron masses and decay 
constantsl^l. On the other hand, earlier attempts to apply xPT predictions for 
the quark-mass dependence of baryon magnetic moments failed, as the higher 
order terms of the chiral expansion quickly dominate the truncated expansion 
as one moves away from the chirai limit. To one meson loop, xPT expresses 
the nucleon magnetic moments asQ 

f^N = f^o + cim^ + C2ml\ogm^ + Csml -\ , (1) 

where /iq and C3 are fitted phenomenologically and ci and C2 are predicted 
by xPT. The TO^logm^ term quickly dominates as to^t moves away from the 
chiral limit, making contact with the lattice results impossible. 

As a resul|L|D£ these early difficulties, lattice QCD results for baryon mag- 
netic momentxH'B remain predominantly based on linear quark mass (or to^) 
extrapolations of the moments expressed in natural magnetons. This approach 
systematically underestimates the measured moments by 10 to 20%. Finite lat- 
tice volume and spacing errors are expected to be some source of systematic 
error. However, xPT clearly indicates the linear extrapolation of the simula- 
tion results is also suspect. It is therefore imperative to find a method which 
can bridge the void between the realm of xPT and lattice simulations. 

We report such a method, which provides predictions for the quark mass 
(or TO^) dependence of nucleon magujetic moments well beyond the chiral limit 
- for more details we refer to Ref. Ej. The method is motivated by studies 
based on the cloudy bag model (CBM), a chiral quark model which preserves 
the correct leading non-analytic behaviour of chiral perturbation theory while 
providing what should be a fairly reliable transition to the regime of large pion 
mass. These studies suggest a Fade approximant which incorporates both the 
leading nonanalytic structure of xPT and Dirac-moment mass dependence in 
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the heavy quark- mass regime. We apply this method to the existing lattice data 
as an illustration of how important such an approach will be to the analysis of 
future lattice QCD calculations. 



2 The Cloudy Bag Model 

The linearized CBM Lagr|aiigian, with pseudoscalar pion-quark coupling (to 
order I/Ztt), is given byO'O 

L = \q{^Y^^, - mq)q - S] 6*^ - ]^qq5s 

111 
+ l^{d^,T:f - -mli:'^ - —q-1^T -TTqSs, (2) 

where B is the bag constant, /tt is the tt decay constant, Oy is a step function 
(unity inside the bag volume and vanishing outside) and 5s is a surface delta 
function. In a lowest order perturbative treatment of the pion field, the quark 
wave function is not affected by the pion field and is simply given by the MIT 
bag model. Our calculation is carried out in the Breit frame with the center- 
of-mass correction for the bag performed via Peierls-Thouless projection. The 
detailed formulas for cajrulating nucleon electromagnetic form factors in the 
CBM are given in Ref . E3. 

In the CBM, a baryon is viewed as a superposition of a bare quark core and 
bag plus meson states. Both the quark core and the meson cloud contribute to 
the baryon magnetic moraents. These two sources are balanced around a bag 
radius, R = 0.7 — 1.1 fmllil. A large bag radius suppresses the contributions 
from the pion cloud, and enhances the contribution from the quark core. The 
CBM reproduces the leading non-analytic behavior of xPT, which has as its 
origin the process shown in Fig. 1(c), along with the other contributions in 
the model. 




B 

(a) (b) (c) 



Figure 1: Schematic illustration of the processes included in the CBM calculation. B and C 
denote intermediate state baryons and include N and A. 

For the ■nNN vertex, we take a phenomenological, monopole form, u{k) = 
{h?~li^)/ (A^-|-fc^), where k is the loop momentum and A is a cut-off parameter. 
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As current lattice simulations indicate thaLjTT,^ is approximately proportional 
to niq over a wide range of quark masses llj, we scale the mass of the quark 
confined in the bag as niq = (m^//i)^ ''m^q \ with m^q^ being the current quark 
mass corresponding to the physical pion mass (/x). In order to obtain a first 
idea of the behaviour within the model between the lowest mass lattice point 
and the experimental data, the parameters i?Oj A and rrSq^ were tuned to 
reproduce the experimental moment while accommodating the lattice data. It 



turns out that the best fits are obtained with 



JO) 



in the range 6 to 7 MeV 



for a bag radius of 0.8fm and A of order 600-700 MeV, which are all quite 
satisfactory. 
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Figvme 2: The proton magnetic moment as calculated in lattice QCD (• LDW Ref. III. ■ WDL 
Ref.B), the cloudy bag model (CBM) and the MIT bag model (MIT). Also illustrated is a fit 
of the simple analytic form given in Eq. (p|) to the CBM results. The point at the physical 
value of is the experimental measurement and is used to constrain the parameters of the 
CBM. 



Inspection of the results of this calculation, in Fig. show clearly that 
the pion cloud contribution to the nuclcon magnetic moments decreases very 
quickly, becoming quite small for large quark masses - especially in the range 
corresponding to the current lattice calculations. In particular, the total pionic 
correction at the first lattice data point (for around 600 MeV) is of the 
order of only 10-15% of the total. Since the present lattice data is based on a 
quenched calculation, which gives incorrect chiral contributions, this is quite 
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good news. It suggests that in using this data as the basis for a correct chiral 
extrapolation to the physical mass quenching should not induce a major error. 

It is also clear from Fig. ^ that below TOtt = 600 MeV the behaviour of the 
magnetic moment is highly non-linear. One clearly needs to account for such 
behaviour in a reliable manner if we are to make believable extrapolations of the 
lattice data. Rather than rely on a model dependent extrapolation method, we 
chose to investigate whether it was possible to find a simple phenomenological 
form, with a sound physical basis, which could do the job. The successful 
conclusion to that search is described next. 

3 Encapsulating Formula 

After considerable effort we found that the following simple Pade approximant, 
which builds in both the linear behaviour in as to^ — > (i.e., a square root 
branch point in fh) and the Dirac moment for large (i.e., ^ I/^T')) was 
able to reproduce the behaviour found in the CBM calculations: 

(0) 

/i^(m^) = -— — „ „ ■ (3) 

1 + -^m^ + (imi 

The fit for the proton case is also shown in Fig. 2. Even more remarkable is 
that the fit parameter a turned out to be quite close to the value required by 
xPT (e.g., a = 4.54 for the proton, compared with a = 4.41 in xPT). 

This result encourages us to propose that Eq.(^, with a taken from xPT, 
should be used as the method of extrapolation in future analyses of lattice data 
for nucleon magnetic moments - and suitably generalized, for all members of 
the baryon octet. As an example, we show in Fig. 3 the result of two parameter 
fits to the proton and neutron data. The nucleon magnetic moments at the 
physical pion mass, obtained from this extrapolation, are 

/ip = 2.85(22) ^lN and ^„ = -1.90(15) /^at, (4) 

which agree surprisingly well with the experimental measurements, 2.793 and 
— 1.913 respectively. We note that the data required to do a fit of the lattice 
results in which covariances are taken into account is no longer available. As 
such, the uncertainties quoted here should be regarded as indicative only. 

4 Conclusion 

We have explored the quark mass dependence of nucleon magnetic moments. 
Quark masses beyond the region appropriate to chiral perturbation theory 
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Figure 3: Extrapolation of lattice QCD magnetic moments for the proton (upper) and 
neutron (lower) to the chiral limit. The curves illustrate a two parameter fit of Eq. (H) to 
the simulation data in which the one-loop corrected chiral coefficient of is taken from 
XPT. The experimentally measured moments are indicated by asterisks. 

have been explored using the cloudy bag model (CBM) which reproduces the 
leading nonanalytic behaviour of xPT while modeling the internal structure 
of the hadron under investigation. We find that the predictions of the CBM 
are succinctly described by a simple formula which reproduces the leading 
nonanalytic behavior of xPT in the limit ttItt —>■ and provides the anticipated 
Dirac moment behavior in the limit ttItt ^ oo. As an example we applied 
this encapsulating formula to the existing lattice data for nucleon magnetic 
moments, leading to surprisingly accurate predictions in comparison with the 
observed values. It will be interesting to see how the fit parameters change as 
finite volume and lattice spacing artifacts are eliminated in future simulations 
and whether the level of agreement seen in this investigation is maintained 
or improved. We strongly advocate the use of the Fade approximant given in 
Eq. (^) in future lattice QCD investigations of octet baryon magnetic moments. 
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